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Abstract. We present a new geometric approach to the 
study of static isolated general relativistic systems for which 
we suggest the name geometrostatics. After describing the 
setup, we introduce localized formulas for the ADM-mass and 
ADM/CMC-center of mass of geometrostatic systems. We 
then explain the pseudo-Newtonian character of these formulas 
and show that they converge to Newtonian mass and center of 
mass in the Newtonian limit, respectively, using Ehlers' frame 
theory. Moreover, we present a novel physical interpretation 
of the level sets of the canonical lapse function and apply it to 
prove uniqueness results. Finally, we suggest a notion of force 
on test particles in geometrostatic space-times. 



1. Introduction 

Static isolated general relativistic systems have been studied from 
a number of perspectives including their regularity, compactification 
and asymptotic considerations, symmetry classifications, construction of 
explicit solutions etc. They serve as models of static stars and black holes. 
Also, they play an important role in R. Bartnik's definition of quasi-local 
mass and his associated conjecture on static metric extensions |12j . 

Here, we present a new geometric approach to the study of static 
isolated systems and their physical properties for which we suggest 
the name geometrostatics. We consider space-times that are static 
(possess a smooth global time-like Killing vector field that is hypersurface- 
orthogonal) and isolated (see below). Static space-times generically possess 
a 3 + 1-decomposition with vanishing shift vector. In this canonical 



decomposition, the canonical lapse function is time-independent and 
coincides with the Lorentzian length of the time-like Killing vector field. 
The spacelike time-slices orthogonal to the time-like Killing vector field are 
all isometric and have vanishing extrinsic curvature, see Figure [I] Their 
induced Riemannian metric is time-independent. We will subsequently 
identify all canonical time-slices. 



Figure 1. The time-slices of a canonically decomposed static space-time. 

For our purposes, static systems are called isolated if the Riemannian 
metric and the lapse function on the time-slice decay suitably fast to the 
flat metric and the constant 1, respectively, at spacelike infinity see [3] for 
a precise definition of our asymptotic flatness condition in the language 
of weighted Sobolev spaces. Moreover, we request that the space-time 
satisfies the vacuum Einstein equations outside some spatially compact 
tube in the space-time (or, in other words, outside some compact set in 
the time-slice). This can be interpreted as a (spatially) finite extension of 
the sources, whether they are matter sources and/or black holes. 

This article is structured as follows: In Section [2] we will introduce the 
central equations of geometrostatics and summarize a few of their central 
analytic properties. In Section[3j we present a novel physical interpretation 
of the level sets of the lapse function of a geometrostatic system and discuss 
some applications of this insight. In Section |4j we will perform a conformal 
transformation into what we suggest to call pseudo-Newtonian variables. 
Moreover, we will define and analyze localized surface integral expressions 
for the mass and center of mass of geometrostatic systems. Finally, in 
Section [5j we will discuss the Newtonian limit of geometrostatics. 

Further details can be found in my thesis [3]. 

2. Geometrostatics 

The Lorentzian metric of a generic static space-time R x M 3 can globally 
be decomposed as 

ds 2 = -N 2 c 2 dt 2 + g, (1) 

where N :— y/—ds 2 (dt, dt) is the (canonical) lapse function arising as the 
Lorentzian length of the time-like Killing vector field dt, c is the speed 



of light, and g is the Riemannian metric induced on the time-slice M 3 . 
Observe that N is non-negative and vanishes only along Killing horizons. 

In the vacuum region outside the matter, the (vacuum) Einstein 
equations imply that these variables satisfy the so-called vacuum static 
metric equations 



where Ric is the Ricci curvature tensor of g, V 2 iV denotes the covariant 
Hessian, and AN denotes the covariant Laplacian of N with respect to g. 
It is well-known that solutions to the vacuum static metric equations are 
real analytic in suitable coordinates [llj . 

We define a geometro static system to be an asymptotically flat 
Riemannian 3-manifold (M 3 ,g) endowed with a smooth positive lapse 
function N so that the vacuum static metric equations ^ and Q are 
satisfied. Hence, geometrostatic systems model the vacuum region outside 
the support of the matter and the horizons of all black holes within a slice 
of an asymptotically flat static space-time. The lapse function N describes 
the lapse of time in the space-time. 

3. The Level Sets of the Lapse 

In Newtonian gravity, the relevant gravitational variable is the Newtonian 
potential. The gradient of the potential defines the force on a unit mass test 
particle. This has a well-known consequence for the equipotential surfaces 
(or level sets of the Newtonian potential): if a test particle is constrained 
into one of these surfaces then the gravitational force does not have a 
tangential component and hence the test particle does not tangentially 
accelerate within the surface, see Figure [5] 



iVRic = V 2 7V 
AN = 0, 



(2) 
(3) 



test body restricted to E 




Figure 2. A test particle constrained to a surface E. 



Surprisingly, the "same" is true for level sets of the lapse function in 
a geometrostatic system. In order to make this rigorous, we make the 
following definitions: Consider a closed smooth surface S C M 3 in a 
geometrostatic system (M 3 ,g,N) arising as the n = no level set of a 
smooth function n : M 3 — s- M. A time-lik^H curve 

M (r) = (t(T),x(T)) (4) 

satisfying x(t) € X is called a tesi particle constrained to E if it is a critical 
point of the time functional 

T(n):= rmr)\+o-(nox(T)-n )}dr, (5) 

J T 

where a g R is a Lagrange multiplier ensuring that all comparison curves 
are also constrained to E. With this notion of constrained test particle, we 
say that a smooth closed surface E is an equipotential surface if every test 
particle constrained to E is a geodesic in E with respect to the induced 
2-metric, see Figure [2j Analyzing the geodesic equation, we find that a 
surface E C M 3 is an equipotential surface in (M 3 ,g,N) if and only if 
E is a level set of N. Thus, the level sets of the lapse function N in 
geometrostatics play precisely the same role as those of the Newtonian 
potential in Newtonian gravity. 

In a static vacuum space-time, the Einstein constraint equations reduce 
to Seal = 0. In particular, the lapse function does not appear in this 
constraint equation. As a consequence, Choquet-Bruhat's theorem (see 
e.g. [3]) on the local existence and uniqueness of solutions to the Einstein 
equations implies that the space-time induced by (M 3 ,g,N) is in fact 
independent of the lapse function N. Combining this view of the lapse 
function with the physical interpretation of its level sets of the lapse as well 
as with the vacuum static metric equations |2]) and ([3| and the assumed 
asymptotic conditions for g and N, we obtain that the lapse function is 
indeed unique if it exists. This is to say that if (M 3 ,g,N) and (M 3 ,g,N) 
are geometrostatic systems, then N — N. We interpret this result as saying 
that "there is only one way of synchronizing time at different locations in 
a geometrostatic space-time such that one sees staticity" just as, for a 
geodesic, "there is only one way of walking along a geodesic such that one 
does not accelerate (up to affine transformations of the curve parameter)". 
The affine freedom of the parameter along the geodesic does not make an 
appearance in the geometrostatic space-time picture because we fixed the 
lapse function to asymptotically converge to 1 at spacelike infinity and 
therewith fixed the time unit. 

1 Here, time-like curves and the time functional are taken with respect to the 
static space-time metric ds 2 = —N 2 c 2 dt 2 + g induced by (M 3 , g, N). 

2 This assumes that the lapse function exists in the first place. 



4. Pseudo-Newtonian Gravity 

The geometrostatic variables g and N are ideal for investigating geometric 
and relativistic effects influencing test particle behavior and the behavior 
of light rays. In order to better understand asymptotic and analytic 
properties of solutions, however, it is more convenient to perform a 
conforms! change and consider the new variables 

7 := N 2 g (6) 
U := c 2 \nN. (7) 

These variables have been used by many authors^] see e.g. [9] . We suggest 
to call them pseudo-Newtonian metric and potential, respectively. The 
vacuum static metric equations ([2]), ^ translate into 

Ric 7 = 2c- A dU®dU (8) 
A 7 f7 = 0, (9) 

where Ric 7 denotes the Ricci curvature tensor of 7 and A 7 denotes the 
7-covariant Laplacian on M 3 . 

N was assumed to converge to 1 asymptotically at spacelike infinity, 
so U must asymptotically tend to 0. Indeed, Kennefick and O'Murchadha 
[5] showed that the asymptotic flatness assumptions incorporated into the 
above definition of a geometrostatic system induce the decay conditions 

7 « = S tJ +0(r- 2 ) (10) 

U = - mC ^+0(r- 2 ) (11) 
r 

as r — > 00 in suitable coordinates. Here, m is the ADM-mass of the slice 
(M 3 ,g), see [TJEIIS]. In [3], we prove asymptotic estimates in weighted 
Sobolev spaces that improve this fall-off result. In particular, we find that 

/ M 2 \ r 2M 2 x l x i , N , , 

7« = (l-^)^ + ^4^+<^ 3 ) (12) 

U = _^_^ll + (r-% (13) 

as r — > 00, where M = mG/c 2 and za £ K 3 is a fixed vector. This 
decay occurs in asymptotically flat 7-harmonic coordinates. As a matter 
of fact, these coordinates coincide with the asymptotically flat (spatial) 
wave-harmonic coordinates on (K x M 3 ,ds 2 = —N 2 c 2 dt 2 +g). The vector 
za can be interpreted as the coordinate vector of the asymptotic center of 
mass of the system, see below. 



3 albeit without explicit reference to the speed of light. 



In terms of the pseudo-Newtonian variables 7 and U and inspired 
by Newtonian gravity, we suggest the following quasi-local definitions of 
pseudo-Newtonian mass and center of mas^ of a geometrostatic system 
(M 3 ,g,N) with associated pseudo-Newtonian variables (7, 17) 

where E is any surface enclosing the support of the matter, v and da are 
the outer unit normal to and area measure of E with respect to 7, and x 
is the vector of 7-harmonic coordinates. 

Surprisingly, both of these expressions are independent of the particular 
choice of surface E (as long as the surface encloses the support of the 
matter, imagine for example a large coordinate sphere). For the pseudo- 



Newtonian mass (14 1, this independence of the surface can be seen by 
combining Equation ^ with the divergence theorem. For the center 
of mass (15), the independence of the surface follows from Equation Q 
combined with Green's formula and the fact that the coordinates are 7- 
harmonic such that A 7 x = 0. We will thus drop the explicit reference to 
the surface when referring to pseudo-Newtonian mass and center of mass. 



Using the asymptotic decay (10), (11), we find that 

m PN = m ADM - (16) 
The total mass of a geometrostatic system is thus localized. It can be read 



off exactly on any surface enclosing the matter. Applying Formula (14) to 
any smooth surface E C M 3 , we immediately obtain a notion of mass for 
an arbitrary part of the system (namely the part bounded by the surface 
E). By the divergence theorem and (pj), the masses of all components in 
a multi-component system add up to the total mass of the system just as 
Newtonian masses do. 



If we combine the asymptotic decay (12), (13) with Formula (15) 
defining the pseudo-Newtonian center of mass, we find that Ipjy = zj±. 
We claim that this vector can indeed be physically interpreted as the 
coordinate vector of the total center of mass of the system. For this, 
we exploit Huang's work [7] showing that the ADM-center of mass pQ 
coincides^] with the CMC-center of mass constructed via a constant mean 
curvature (CMC) foliation near infinity by Huisken and Yau [5] and 



4 We note that our discussion of center of mass only applies to systems with 
non-vanishing mass. 

5 under precise fall-off conditions at spacelike infinity that are satisfied here. 



generalized by Metzger [10]. Using again (12 1 and (13), we obtain 



ZPN = ZA = ZADM = ZCMC 



(17) 



which justifies the name center of mass for the quantities zpn and za ■ 
Moreover, this shows that the center of mass of a geometrostatic system 
is also localized. As above, we obtain a notion of center of mass for an 
arbitrary part of the system (namely the part bounded by the surface S) . 
By Green's formula, 7-harmonicity of the coordinates, and ([9]), the centers 
of mass of all components in a multi-component system add up to the total 
center of mass of the system just as Newtonian centers of mass do. 



5. The Newtonian Limit of Geometrostatics 

Intuitively, the mass and center of mass of a relativistic system should 
converge to the Newtonian mass and center of mass of its Newtonian limit 
c — > 00. To make this precise, we use Ehlers' frame theory [6j which unifies 
general relativity (GR) and Newton-Cartan gravity (NC) into a common 
geometric framework with geometric variables g, h, T and matter tensor T. 

In frame theory, taking the Newtonian limit corresponds to taking 
a parametric curve of solutions of GR with parameter A = c -2 to its 
limit A —> 0, see Figure [3j Modeling Killing vector fields, staticity, and 
asymptotic flatness in frame theory, we show that the pseduo-Newtonian 
potential converges to the Newtonian potential and the metric 7 converges 
to the flat metric along any family of geometrostatic systems that possesses 



a static Newtonian limit. As the localized pseudo-Newtonian formulas ( 14 ) 



and ( 15 ) are nearly identical with the Newtonian formulas, this proves that 



indeed the relativistic mass and center of mass converge to their Newtonian 
counterparts. 



g,h,T,T 




g-.h.T.T /K 




Figure 3. The universe of Ehlers' frame theory and the Newtonian limit. 



6. Conclusion and Outlook 

Geometrostatic systems share many features with Newtonian ones. First, 
the level sets of the lapse function N (or, equivalently, those of the pseudo- 
Newtonian potential U — c 2 In N) have the same equipotential properties 
as the level sets of the Newtonian potential. We thus define the force on a 
unit mass test particle as F :~ — V 7 Z7 where V 7 denotes the 7-covariant 
gradient. A second Newtonian type law holds for this notion of force [3]. 

Secondly, the total mass and center of mass of a geometrostatic system 
are localized. We put forward explicit geometric formulas for them that 
also allow for the computation of the notions of mass and center of mass 
of individual regions. We applied these formulas to prove consistence of 
ADM-mass and ADM/CMC-center of mass with the Newtonian limit. 

This fact and the notion of force might turn out useful for the study of 
the well-known static n-body problem and of Bartnik's conjecture |12j . 
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